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ABSTRACT 

We present the first measurement of the relationship between the Sunyaev-Zel'dovich effect (SZE) 
signal and the mass of galaxy clusters that uses gravitational lensing to measure cluster mass, based 
on 14 X-ray luminous clusters at z ~ 0.2 from the Local Cluster Substructure Survey. We measure 
the integrated Compton y-parameter, Y, and total projected mass of the clusters (A/ql) within a 
projected clustercentric radius of 350 kpc, corresponding to mean overdensities of 4000—8000 relative 
to the critical density. We find self-similar scaling between Mgl and Y, with a scatter in mass at 
fixed Y of 32%. This scatter exceeds that predicted from numerical cluster simulations, however, it is 
smaller than comparable measurements of the scatter in mass at fixed Tx- We also find no evidence 
of segregation in Y between disturbed and undisturbed cl usters, as had been seen with Tx on the 
same physical scales. We compare our scaling relation to the iBonamente et al.1 relation based on mass 
measurements that assume hydrostatic equilibrium, finding no evidence for a hydrostatic mass bias in 
cluster cores (Mql = 0.98 ±0.13 Mhse), consistent with both predictions from numerical simulations 
and lensing/X-ray-based measurements of mass-observable scaling relations at larger radii. Overall 
our results suggest that the SZE may be less sensitive than X-ray observations to the details of cluster 
physics in cluster cores. 

Subject headings: cosmology: observations — galaxies: clusters: general — gravitational lensing 
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1. INTRODUCTION 

The Sunyaev-Zel'dovich effect (SZE) is a weak distor- 
tion of the cosmic microwave background (CMB) spec- 
trum introduced as CMB photons propagate through 
foreg round galaxy clusters ijSunvaev &: Zel'dovichl 11970 . 
1972). Because this signal does not suffer cosmological 
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dimming and is expected to closely track cluster mass, 
the SZE is a potentially powerful tool for producing large, 
mass-limited galaxy cluster samples that can be used to 
constrain dark energy, under the proviso that the rela- 
tionship betwee n cluster mass and the observed signal is 
well calibrated (|Carlstrom et al.ll2002l ). 

The relationship between the SZE observable, the in- 
tegrated Compton y-parameter, Y, and cluster mass is 
difficult to calibrate because of the difficulty of mea- 
suring mass directly. To date all examinations of the 
mass-V relationship have assumed hydrostatic equilib- 
rium (HSE) when converting SZE or X-ray measure- 
ments of the intr acluster gas into estimates of the total 
mass (Mhse; e.g. iMorandi et al1l2007t IBonamente et al.1 
I2008h . These methods measure scaling relations between 
closely related quantities derived from the same data; 
deviations in one parameter will therefore inevitably be 
correlated with the other, and the true intrinsic scat- 
ter will likely b e underestimated. Furthermore, both 
simulations (e.g. iNagai et al. [20071 iJeltema et ail 120081 : 
iPiffaretti fcVa ldarnini 2008F Lau et al.ll2009h and obser- 
vations (jMahdavi et al.ll2008h suggest that the HSE as- 
sumption is incorrect; non-thermal sources of pressure 
support may bias A/hse or increase the scatter with Y. 

The scatter in mass at fixed Y is an important quan- 
tity to determine in preparation for the SZE cluster sur- 
veys designed to constrain the dark energy equation of 
state, as scatter in the mass-observable relation is gener- 
ally degen erate with the cosmological paramete rs of in- 
terest fe.g. lSmith eVal1l2CMISt"anek et al.ll2006h . These 
surveys will rely on self-calibration techniques to build 
cluster mass functions from SZE-selected cluster cata- 
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logs, however these techniques can be rendered ineffec- 
tive if the scatter is mu ch larger than the few percent 
predicted by simul ations (|da Silva et al.ll2004l ; lMotl et al.1 
2005; Nagai 2006). At present there are only weak ob- 
servational limits, but independent measurements of the 
mass-F scatter ca n be used to refine the self-calibration 
(|Lima fc Hull2005Tl . 

The Local Cluster Substructure Survey (LoCuSS 16 ; G. 
P. Smith et al. 2009, in preparation) is a morphologi- 
cally unbiased survey of X-ray luminous galaxy clusters 
in a narrow redshift range (0.15 < z < 0.3). The sur- 
vey unites a variety of probes of total mass, intracluster 
gas, and cluster member galaxies (e . g..lSmith et alJl2005L 
hereafter [H; IZhang et all 120081: lOkabe et alj l2009f ). 
One of the main goals of the survey is to measure 
the shape, normalization and scatter of cluster mass- 
observable relations as an input to cosmological cluster 
studi es. Here we use a pilot s a mple of clusters, assembled 
from IS051 . iBonamente et alj (|2006l ). and early LoCuSS 
observations with the Sunyaev-Zel'dovich Array (SZA), 
to compare the SZE signal with cluster masses derived 
from gravitational lcnsing (Mql) and thus to construct 
the first mass-Y relation that is independent of the as- 
sumption of HSE. We summarize the observations and 
modeling in Section [2] The Mq^—Y relationship, scat- 
ter, and its implications are discussed in Section [3J We 
assume Q.m—0-3, 0.^=0.7, and h=0.7; the median clus- 
ter redshift is z — 0.222, at which 1" corresponds to a 
physical scale of 3.6 kpc. 

2. OBSERVATIONS AND MODELING 

We analyze SZE and gravitational lensing observa- 
tions of 14 clusters, listed in Table U SZE results 
for eight of these come from 



IBonamente et all (2 006. 200J )_, the observational details 

New SZE measurements 



being found in lLaRoque et a 



LaRoaue et alj (|2006h and 



of A 209, A 383, A 1763, A 2219, A 2537, and Z 2701 were 
obtained at 31 GHz using the SZA, which has baselines 
of 350-1300A for SZE sensitivity (angular scales of 110"- 
7') and 2-7.5 kA for radio source removal. A detailed de- 
scription of SZA observations and analysis methods can 
be found in lMuchovei et all (|2007f h Typical integration 
times were ^ 30 — 40 hr per cluster, and the rms noise 
per beam ~ 0.2mJy. In six of 14 clusters we detect radio 
sources within an arcmin ute of the cluster c enter, these 
sources were presented in lCoble et al.1 (|2007l ). 

Our interferometric observations do not measure di- 
rectly the total SZE flux of the cluster within an aper- 
ture, this must be derived from a model cluster pro- 
file. As in previous works, we have modeled these 
clusters with an isothermal /3-model, which has been 
shown to provide Y measurements indistinguishable from 
those from a simulation -motivated pressure prof ile for 
radii smaller than r V i r /3 (|Mroczkowski et aLl l20091 . The 
shape parameters for the /3-model and an isothermal X- 
ray temperature were jointly derived from the SZE data 
and from Chandra X-ray image s after applying a 100 kpc 
core cut (|LaRoque et al.l [20061 ) . Best-fit cluster parame- 
ters were obtained using the Markov chain M onte Carlo 
method described in lBonamente et alj (|2004f ). The fluxes 
of detected radio sources a nd central sources p resent in 
the 1.4 GHz NVSS catalog ([Condon et al.|[l998h but un- 

16 http:/ /www. sr. bham.ac.uk/locuss 



detected at 31 GHz are included as parameters in the 
Markov chains and marginalized over. 

Details of the gr avita tional lensing mass measure ments 
are described by IS05I and iRichard et alj |2007); the 
source of each measurement is listed in Table [TJ Briefly, 
the mass measurements are based on Hubble Space Tele- 
scope (HST) imaging of the cluster cores, ground-based 
spectroscopy of multiply imaged galaxies identified in 
these data, and by weakly sheared background galax- 
ies. These data constrain parameterized models of the 
projected mass distribution in the cluster cores compris- 
ing one or more cluster-scale mass components (to de- 
scribe the dark matter and intracluster gas) plus typ- 
ically ~ 30 galaxy-scale mass components per cluster. 
Measurements of the projected cluster mass and the un- 
certainty on the mass are then obtained by integrating 
the mass distributions of the family of models within the 
chosen confidence interval in parameter space (in this let- 
ter we quote errors at 68% confidence) out to the chosen 
r adius . 

IS05l (see also lSmith fc Tavlor|[200l ) used the structure 
of the cluster mass distributions inferred from these lens 
models - the substructure fraction, / su b - in conjunction 
with Chandra observations to classify the cluster cores 
as "disturbed" or "undisturbed". The most straightfor- 
ward criterion in this classification was the offset between 
the peaks of the X-ray emission and the lensing-based 
mass map, in the sense that a statistically significant off- 
set indicates that the cluster core is disturbed in some 
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2008) 



We adopt the classification of 



five additional clusters in our sample (Table Q]). 

3. RESULTS AND DISCUSSION 

3.1. The Mass-Y Scaling Relation 

We first define the aperture within which A/gl and 
Y are measured. After some e xper imentation we chose 
to retain the aperture used by IS05I . a clustercentric ra- 
dius of 250h~ 1 kpc, or 350 kpc in our adopted cosmol- 
ogy. This fixed physical aperture is well-matched to the 
2JST/WFPC2 and Advanced Camera for Surveys (ACS) 
fields of view, and is somewhat larger than the typical 
resolution of our SZE data. The shape of our /3-model fit 
is jointly determined from both SZA and high-resolution 
X-ray data and thus is well resolved. Mql and Y values 
measured at this radius are listed in Table Q] and plotted 
in Figure [T] 

We fit our projected mass (Mql) and Y data points in 
the log(y£>^)-log(MGL) plane, where a power law scal- 
ing relation takes the form log (Mql) — ct+f3 log (YD\) . 
Starting with the self-similar scaling (Kaiser 1986), it can 
be shown that the relation between mass and Y within a 
fixed physical radius has a slope = 0.5 with no depen- 
dence on the redshift evolution of the Hubble parameter. 
We seek to measure the slope and normalization of the 
scaling relation, as well as the intrinsic scatter between 
these parame t ers, a nd therefore follow the prescription of 
IWeiner et all (|2006l ) for linear regression in the presence 
of intrinsic scatter. This Baycsian method includes the 
intrinsic scatter in the y-coordinate (a y ) as a parameter 
of the fit and maximizes the probability of the model (a, 
(3, <jy) given the data and errors. The likelihood takes 
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TABLE 1 

Cluster Sample 



Cluster 


Kcasmit 


Y 


v n2 


b/jLj 




Lensmg 


Classification 






(10 1U ) 


(10~ 5 Mpc 2 ) 


Rcf. 


(1O 14 M ) 


Ref. 




A 68 


0.255 




Q fS7-(-n ^9 






2 3 


TllQ't"11"V"r~\f~>/H 
J_vl£> L HI UCU 


A 209 


0.206 


n n ,1-1-0 1 A 


4 ^7-1-0 fi7 


4 




2 


TllQ't~11"V"r~\f~>/H 


A 267 


0.230 


n c;q+o Ofi 


q 08+0 34 




2 20+0 '34 


j 


Til cM 1 CnPn 
i_/ lo U 111 k_/^v_l 


A 383 


0.188 


0.39±0.04 


1.61±0.16 


4 


3.71±0.82 


2 


Undisturbed 


A611 


0.288 


0.39±0.04 


3.13±0.34 


1 


2.12±0.05 


5 


Undisturbed 


A 773 


0.217 


1.03±0.11 


5.40±0.57 


1 


4.03±0.12 


2.5 


Disturbed 


Z2701 


0.214 


0.28±0.03 


1.46±0.16 


4 


1.92±0.07 


5 


Disturbed 


A 1413 


0.143 


1.83±0.26 


4.90±0.70 


1 


2.59±0.50 


5 


Undisturbed 


A 1689 


0.181 


1.86±0.15 


7.51±0.60 


1 


7.44±0.05 


6 


Disturbed 


A 1763 


0.288 


0.56±0.06 


3.10±0.32 


1 


1.42±0.54 


2 


Disturbed 


A 1835 


0.253 


1.03±0.07 


6.82±0.48 


1 


3.35±0.06 


2.5 


Undisturbed 


A 2218 


0.171 


1.12±0.10 


4.23±0.38 


1 


4.23±0.09 


2 


Disturbed 


A 2219 


0.228 


1.12±0.05 


6.27±0.26 


4 


3.48±0.07 


2 


Disturbed 


A 2537 


0.297 


0.42±0.03 


3.47±0.24 


4 


1.75±0.90 


7 


Disturbed 



References. — (1) Bonamente ct al. (2008); (2) S05; (3) Richard ct al. (2007); (4) this work; (5) Richard et al. 
"all H2007I) ; 



2009, in preparation; 

(6) [Limousin ct al. (2007]); (7) Hamilton-Morris ct al. 2009, in preparation. 

Note. — Y and Mql measured at 350 kpc radius. Morphological classification from lensing and X-ray data. See S05 and Smith 
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Fig. 1. — Scaling between projected mass (Mql) based on grav- 
itational lens models and Y, both measured within a projected 
clustercentric radius of 350 kpc. Clusters classified as undisturbed 
in Table [T] are shown as filled circles, disturbed clusters are shown 
as open boxes. The solid line shows the scaling relation when the 
slope is fixed to the self-similar value (/3 = 0.5), the dashed line is 
obtained when the slope is also a free parameter. 



the form (Equation A3 of IWeiner et al.i r2006) 

\Vi - (a + f3xi)) 2 



(1) 



where the ei are the x or y errors in the individual data 
points. 

We perform the regression in two different ways: first, 
with both the normalization (a) and slope (/?) of the 

fit free we obtain a = W.btl l and P = °- 47 -o.20> 
with an intrinsic scatter in mass at fixed Y of <Jm\y = 
32 ± 4%. With the slope fixed to the self-similar value 
{(3 = 0.5) we obtain a = 16.68 ± 0.04, again with 
a M\Y = 32 ± 4%. Despite large uncertainties, the 
Mgl — y sl ope is consistent with self-s imilarity. For 
comparison, iBonamente et~a l. (2008) and Mo randi et ail 



(|2007f ) found that the A/hse — Y scaling relation is con- 
sistent with self-similarity at ^2500 (the radius where the 
average interior density is A=2500 times the critical den- 
sity of the universe, typically 30—100% larger than the 
radius used here). In contrast. IS05I were unable to con- 
strain the slope of the Mq^—Tx relation, suggesting that 
the intrinsic scatter in Mq^-Y is smaller than in Mql- 
Tx- 

3.2. Intrinsic Scatter 

We now compare the intrinsic scat ter in the Mq^—Y 
relation with that in Mgl, - Tx from S05]. We first use 
our regression method to re-fit their data, taking the op- 
portunity to adopt the theoretical slope (Tx = BMq^ , 
f3 = 1, in their nomenclature) appropriate for this scaling 
relation within a fixed aperture. With their 10 clusters 
we find a scatter in mass of <Tm\t = 41±5%, which is 28% 
larger (1.5a significance) than <Jm\y- Indeed, both (J m\y 
and <Jm\t are large compared to the scatter predicted 
from numerical simulations (e.g. Nagai 2006), likely due 
to a combination of modeling uncertainties, astrophysi- 
cal processes in cluster cores, and projection effects. The 
effects of modeling uncertainties are hard to quantify but 
should not be ignored, as both Mql and Y are derived 
from parameterized models. To the extent that these 
models inadequately represent the cluster or underes- 
timate the uncertainty in the derived parameters, our 
inferred scatter will be artificially increased. Thus our 
measured Cm|y is an upper limit to the true scatter be- 
tween mass and Y in cluster cores. 

In general, cluster mass measurements based on grav- 
itational lensing and SZE observations are more suscep- 
tible to projection effects than those based on X-ray 
data because the density-squared dependence of the X- 
ray emissivity limits contributions from material along 
the line of sight. Never theless, simulations (jHolder et al.1 
l2007t IShaw et aill2008h suggest that the contribution to 
the SZE signal from projection is unimportant for clus- 
ters with M V i r > 1O 14 /i _1 M0, with the projection effects 
diminishing further at the small radii used here. Grav- 
itational lensing is sensitive to all mass along the line 
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of sight through the cluster, which will likely increase 
the scatter in our scaling relation, however the contribu- 
tion of correlated large-scale structure to lensing-based 
mass measurements should decrease at smaller radii in a 
manner similar to the impact on the SZE signal. We also 
emphasize that our sample is X-ray selected with no con- 
sideration paid to the presence/ absence of strong lensing 
arcs, so it should not suffer the halo orientation bias that 
boos ts the core mass in s trong-lensing-selected samples 
(e.g. iHennawi et al.ll2007t) . Indeed, 5/14 clusters con- 
tain no obvious strong-lcnsing signal. Projection-induced 
scatter therefore appears insignificant for Y and may be 
somewhat more important for Mql, although this affects 
both <7m\y and "MIT- 
Given that the mass measurement methods and 9/14 
of the clusters are identical between this Letter and IS051 
we expect the difference between the scatter in these 
two relations to be dominated by the relative sensitiv- 
ity of Y and Tx to the gas physics of cluster cores. 
The intracluster medium (ICM) in cluster cores is of- 
ten dist urbed by active ga lactic nucleus (AGN) activ- 
ity (e.g. iFabian et al.ll2006f). an d cluster-cluster mergers 
(jMarkevitch &; Vikhlininl 120071) . For this reason, clus- 
ter cores are often excluded when de riving X-ray qbserv- 
ables for use as mas s proxies (e.g. lArnaud et alj 120051 ; 
iKravtsov et"afll2006t iMaughanl 120071). We therefore as- 
cribe the difference between <Jm\y and om\t to the dif- 
ference in the SZE and X-ray emissivity dependencies on 
ICM density and temperature. The SZE is unaffected 
by isobaric disturbances (excluding relativistic compo- 
nents such as cosmic rays), while the energy-integrated 
X-ray emissivity varies approximately as the square of 
the de nsity, which may be sig nificantly perturbed in the 
core. iPfrommer et al.l (|2007l ) found potentially impor- 
tant changes in the SZE and X-ray signals from cool cores 
in simulations incorporating cosmic rays, but the effect 
and therefore the induced scatter was found to be larger 
for X-ray emission. In larger measurement apertures, 
such as the frequently used r25oo and rsoo (350 kpc corre- 
sponds to A=4000— 8000 for our sample), core-removed 
X-ray mass pro xies are found to corr elate more tightly 
with mass (e.g. IVikhlinin et al.ll2009f ). We expect that 
°~m\y "will also decrease when measured on these scales. 

3.3. Structural Segregation 

IS05I found that the normalization of the Mq-^—Tx re- 
lation for disturbed clusters was ~ 40% hotter than for 
undisturbed clusters at ~ 2.5<r (after excising the central 
region of the X-ray data). We re-fit the Mq^-Y relation 
to the undisturbed and disturbed sub-samples (Table [T]) 
finding no significant difference in normalization be- 
tween undisturbed and disturbed clusters (a U ndisturbed — 
o-disturbed = 0.0±0.1, (3 fixed to 0.5). Again, this suggests 
that the integrated SZE signal, even in cluster cores, is 
less sensitive to the ICM physics than the X-ray temper- 
ature, although the significance of t his re sult is low given 
the small samples used here and in lS05l 

3.4. Hydrostatic Equilibrium 

As discussed in Sections HI and I3~2l the ICM may be far 
from HSE; we search for this e ffect in our data by com - 
bining our results with those of lBonamente et"afl (|2008l ) 
who studied the Mhse-^ relation. Our lensing measure- 
ments are sensitive to the projected mass (M cy i); we 



therefore convert them to the spherical masses (M sp h) 
required for comparison with the hydrostatic analysis. 
We do this by assuming that the halo mass distribu- 
tions ar e described by the N avarro- Frenk- White (NFW) 
profile ([Navarro et al.l [l996). which is specified by mass 
and concentration, M V i r and c V i r ; given Mql and c V i r 
we can than derive an approximate conversion factor 
M cy i/M sp h at 350 kpc. This conversion factor depends 
sensitively on c V i r , but there is considerable theoretical 
uncertainty about the mean value and scatter in con- 
centration at fixed mass, and its variation with redshift 
and cluster mass. Without direct measurements of c V i r 
for most of our cluste rs we margin a lize o ver the range 
of values observed by lOkabe et al.l (|2009f ) for 19 clus- 
ters from the LoCuSS sample, which have very similar 
masses and redshifts to our samp le. We randomly ch oose 
a concentration from the 19 in lOkabe et alj (|2009t ) for 
each of our clusters, calculate the corresponding value 
of MGL,s P h=M G i J xM S ph/M C yi, and fit a scaling rela- 
tion between the estimated McL.sph and Y. We re- 
peat this process thousands of times, and derive a mean 
normalization of the McL.sph— Y relation of aGL.spfc = 
16.42 ± 0.05 for j3 = 0.5. This value can be directly com- 
pared to the normalizati on of the Mhse - Y relation from 
iBonamente et al.l (|2008l) : a H SE = 16.425 ± 0.016. The 
statistically insignificant difference between these values 
implies a mass ratio of MGL,s P h/-^HSE=0.98 ± 0.13. 

Using hydrostatic masses from X-ray obs ervations and 
weak lensin g masses from t he lit erature, IZhang et al.l 
(|2008| ) and IVikhlinin et al.l (|2009f ) estimate the weak 
lensing to hydrostatic mass ratio to be 1.09 ± 0.08 
and 1.01 ± - 11, respectively, measured at ^00 ■ 
iMahdavi et al.l (|2008f ) . however, demonstrate a radial 
trend in this ratio between r25oo and rsoo, suggestive of 
a 20% deficit in hydrostatic masses at rsoo that cannot 
be accounted for by the projection of unrelated struc- 
tures along the line of sight. Our result, interior t o 
the innermost radius examined by IMahdavi et a l. (2008), 
is consistent with their non-detection of a bias in the 
hydrostatic mass at small radii. The precision of our 
measurement is limited by the intrinsic scatter in our 
relation, with a smaller contribution coming from the 
unknown halo concentration values required to convert 
between Mql and McL,sph- Future combinations of 
weak and strong lensing measurements, specifically joint 
strong+weak lens modeling of the clusters combining the 
Subaru and HST data (G. P. Smith et al. 2009, in prepa- 
ration) should allow us to evaluate the bias as a function 
of radius. 

4. SUMMARY 

We have presented the first calibration of the mass-F 
relation based on gravitational-lensing measurements of 
cluster mass, based largely on previously published SZE 
data and gravitational lens models. We construct the 
relation at a fixed physical radius of 350kpc to minimize 
uncertainties arising from extrapolation of cluster mod- 
els based on both datasets. In cont rast to the mass-Tx 
relation within the same aperture (|S05l ). we succeed in 
fitting a model with both free slope and normalization. 
The best-fit slope is consistent with the self-similar pre- 
diction, and the intrinsic scatter in mass at fixed Y is 
32%, in contrast to the 41% scatter in mass at fixed Tx- 
We also fit the relation to sub-samples of disturbed and 
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undisturbed cluster cores but find that the best-fit nor- 
malizations for these sub-samples are consistent within 
the errors. Finally, we c ombine our results with those of 
iBonamente et al. (2008) to test whether the cluster cores 
are in HSE. Cluster core masses estimated from lensing 
and SZE data (assuming HSE) are, on average, consis- 
tent within the errors, suggesting that departures from 
equilibrium arc modest in cluster cores. We conclude 
that the difference between the mass-l" and mass-Ty re- 
lations is mainly attributable to the relative insensitivity 
of the SZ effect to the physics of the ICM in cluster cores. 
Future articles will explore the mass-F relation at larger 
radii through weak lensing and in larger samples. 
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